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Abstract 

We review several results related to the problem of a quantum particle in a random 
environment. 

In an introductory part, we recall how several functionals of the Brownian motion arise 
in the study of electronic transport in weakly disordered metals (weak localization). 

Two aspects of the physics of the one-dimensional strong localization are reviewed : some 
properties of the scattering by a random potential (time delay distribution) and a study of 
the spectrum of a random potential on a bounded domain (the extreme value statistics of 
the eigenvalues). 

Then we mention several results concerning the diffusion on graphs, and more generally 
the spectral properties of the Schrodinger operator on graphs. The interest of spectral 
determinants as generating functions characterizing the diffusion on graphs is illustrated. 

Finally, we consider a two-dimensional model of a charged particle coupled to the random 
magnetic field due to magnetic vortices. We recall the connection between spectral properties 
of this model and winding functionals of the planar Brownian motion. 

PACS numbers : 72.15.Rn ; 73.20.Fz ; 02.50.-r ; 05.40. Jc. 
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1 Introduction 

1.1 Weak and strong localization 

At low temperature, the electric conductivity a of metals and weakly disordered semiconductors 
is determined by the scattering of electrons on impurities. It is given by the Drude formula 

n e e 2 T e 

oo = , 1 

m 

where e and m are the charge and mass of the electron respectively. n e is the electronic density 
and r e the elastic scattering time 1 . This purely classical formula is only valid in a regime where 
quantum mechanical effects can be neglected. This is the case if the elastic mean free path 
of electrons t e = vpr e is large compared with the De Broglie wave length Xp = 2Trh/mvp 
corresponding to the Fermi energy (vp is the Fermi velocity). 

Strong localization. - When these two length scales are of the same order £ e ~ Xp (strong disorder) 
the fact that the electrons are quantum objects must be taken into account and the wave like 
character of these particles is of primary importance. It is indeed this wave character which is 
responsible for the localization phenomenon. The multiple scattering on impurities distributed 
randomly in space creates random phases between these different waves which can interfere 
destructively. These interference effects reduce the electronic conductivity. In the extreme case 
of very strong disorder the waves do not propagate anymore and the system becomes insulating. 
This is the strong localization phenomenon which was conjectured by Anderson in 1958. 

Weak localization. - In the 80's it was realized that even far from the strong localization regime 
the quantum transport is affected by the disorder. Diagrammatic techniques used in the weak 
disorder limit Ap <C £ e , were initiated by the works of AFtshuler, Aronov, Gor'kov, Khmel'nitzkii, 
Larkin and Lee [1161 112j (see ref. ^3] f° r an introduction and ref. jl] for a recent presentation) . 
In this regime, called the weak localization regime, the Drude conductivity gets a small sample 
dependent correction whose average, denoted (Act), is called the "weak localization correction". 
(• • •) denotes averaging with respect to the random potential. From the experimental side, this 
phenomenon is well established and has been the subject of many studies (see refs. |.'S4( llfij for 
review articles). 

1 The elastic scattering time r e is the time characterizing the relaxation of the direction of the momentum of the 
electron. The conductivity is proportional to r e for isotropic scattering by impurities. For anisotropic scattering 
the conductivity involves a different time r tr called the "transport" time I17| (see ref. |1] for a discussion within 
the perturbative approach). 



2 



Phase coherence and dimensional reduction.- The localization phenomenon comes from the 
interplay between the wave nature of electronic transport and the disorder. This manifestation 
of quantum interferences requires that the phase of the electronic wave is well defined, however 
several mechanisms limit the phase coherence of electrons in metals, among which are the effect 
of the vibrations of the crystal (electron-phonon interaction) or the electron-electron interaction. 
We introduce a length scale L v , the phase coherence length, that characterizes the length over 
which phase breaking phenomenon becomes effective. The lack of phase coherence in real systems 
is the reason why the strong localization regime has not been observed in experiments on metals. 
In dimension d = 3 the strong localization regime is only expected to occur for sufficiently strong 
disorder 2 . However, strong localization can also be observed in the weak disorder limit (Xp <C £ e ) 
by reducing the dimensionnality 3 . It has been shown in the framework of random matrix theory 
that the localization length of a weakly disordered quasi-ld wire behaves as A ~ N c £ e , where N c 
is the number of conducting channels 4 |791I8U1 11661 13U| . Therefore the strong localization regime 
is expected to occur when coherence is kept at least over a scale A. At low temperature, in the 
absence of magnetic impurity, phase breaking mechanisms are dominated by electron-electron 
interaction, which leads to a divergence of the phase coherence length L V (T) oc (A^/J 1 ) 1 / 3 
predicted in ref. [3] and verified in several experiments like |224U213"ll84| (see ref. |180| for recent 
measurements down to 40mK). Therefore the temperature below which strong localization might 
be observable is given by L V (T*) ~ A, which leads to T* ~ l/(N^dr e ). The crossover temperature 
in metals is out of the experimental range, however it becomes reachable in wires etched at the 
interface of two semiconducting materials, when the number of conducting channels is highly 
reduced, and the manifestation of the strong localization has been observed in ref. |1U9| . 

Strictly one- dimensional case. - In a weakly disordered and coherent quasi-ld wire, the weak 
localization regime only occurs for length scales intermediate between the elastic mean free path 
and the localization length £ e <C L -C A. In the strictly one dimensional case, since A ~ 4£ e 
for weak disorder 5 , such a regime does not exist and the system is either ballistic (£ e S> L) or 
strongly localized (£ e <C L). 

Anderson localization in one dimension has been studied by mathematicians and mathematical 
physicists from the view point of spectral analysis and in connection with limit theorems for 
products of random matrices [IS]. A breakthrough was the proof of wave localization in one 
dimension by Gol'dshtein, Molchanov & Pastur |113j . Although some progress has been made, 
the multidimensional case is still out of reach and the subject of weak localization has almost 
not been touched in the mathematical literature. One of the main fields of interest in the last 
twenty years is the investigation of random Schrodinger operators in the presence of magnetic 
fields (see for example the refs. |22U1 1461 1771 1781 11UU1 HUlj ; for recent results on Lifshitz tails 
with magnetic field see refs. |lfll[ 1153] ). From the physics side significant advances have been 
realized. Field theoretical methods based on supersymmetry provide a general framework for 

2 Following the scaling ideas initiated by Thouless |156| and Wegner |221l I219| it was shown in ref. that 
the localization-delocalization transition exists only for dimension d 3. In d = 1 and d — 2, the fully coherent 
system is always strongly localized, whatever the strength of the disorder is. 

3 The effective dimension of the system is obtained by comparing the sample size with the phase coherence 
length L v . For example, a long wire of length L and of transverse dimension W is effectively in a Id regime if 
W < L v < L. 

4 The localization length predicted by the random matrix theory (RMT) can be found in ref. |3U| : Armt = 
[/3(N C — 1) + 2]£ e , where /3 — 1, 2, 4 is the Dyson index describing orthogonal, unitary and symplectic ensembles, 
respectively. Note that one must be careful with the coefficient involved in this relation since the definitions of 
Armt and £ e differ slightly in RMT and in perturbation theory. 

5 The elastic mean free path £ e = vft b is given by the self energy l/(2r e ) = -ImE R (£). For example, for a 
Gaussian disorder with local correlations, (V(x)V(x')) = wS(x — x'), we obtain l/r e ~ 2npow for a weak disorder, 
where po is the free density of states. In one dimension po = 1/(tvvf), therefore £ e ~ v%/(2w), which coincides 
with the high energy (weak disorder) expansion of the localization length A ~ 2v F /w 1 151 1157| . 
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disordered systems and also allow establishing some links with quantum chaos 86 . 

While writing this review, we have tried to collect a large list of references which is however 
far from being exhaustive. Ref. |157j provides a reference book on strong localization, mostly 
focused on spectrum and localization properties (see also ref. |158j for a review on Id discrete 
models and Lifshitz tails). A recent text about disorder and random matrix theory is ref. |86| . 
Many excellent reviews have been written on weak localization among which |34| I55j (for the 
role of disorder and electron-electron interaction, see the book |87| or |152j ). A recent reference 

is en 

1.2 Overview of the paper 

Section |^1 shows that several functionals of the Brownian motion arise in the study of elec- 
tronic transport in weakly disordered metals or semiconductors (weakly localized). The brief 
presentation of weak localization given in sections 12 . 1 1 and 12 . 21 follows the heuristic discussion of 
refs. [131 155j . which is based on the picture proposed by Khmeknitzkh & Larkin. In spite of its 
heuristic character, it allows drawing suggestive connections with well known functionals of the 
Brownian motion. 

The sections 01 and 0] deal with problems of strong localization in one dimension. A pow- 
erful approach to handle such problems is the phase formalism 6 (presented in refs. [151 1157j 
for instance). This formalism leads to a broad variety of stochastic processes. Section El dis- 
cusses scattering properties of a random potential and section |I] studies spectral properties of a 
Schrodinger operator defined on a finite interval. 

In section we review some results obtained for networks of wires (graphs), which can be 
viewed as systems of intermediate dimension between one and two. We will put the emphasis 
on spectral determinants, which appear to be an efficient tool to construct several generating 
functions characterizing the diffusion on graphs (or its discrete version, the random walk). 

Finally, in section H3 we show that the physics of a two dimensional quantum particle sub- 
mitted to the magnetic field of an assembly of randomly distributed magnetic vortices involves 
fine properties of the planar Brownian motion. 

2 Functionals of the Brownian motion in the context of weak 
localization 

2.1 Feynman paths, Brownian motion and weak localization 

In the path integral formulation of quantum mechanics, each trajectory is weighted with a phase 
factor e lS l h where S is the classical action evaluated along this trajectory. The superposition 
principle states that the amplitude of propagation of a particle between two different points is 
given by the sum of amplitudes over all paths connecting these two points. This formulation 
of quantum mechanics is most useful when there is a small parameter with respect to which 
one can make a quasiclassical expansion. One encounters a similar situation in the derivation 
of geometrical optics starting from wave optics. In this case the small parameter is the ratio 
of the wave length to the typical distances which are involved in the problem and the classical 
paths are light rays. In the context of weak localization the small parameter is Ai?/-£ e and the 
quasiclassical approximation amounts to sum over a certain subset of Brownian paths. The fact 
that Brownian paths come into the problem is not so surprising if one goes back to our previous 
physical picture of electrons performing random walk due to the scattering by the impurities. 

6 The phase formalism is a continuous version of the Dyson- Schmidt method I83II188| . A nice presentation can 
be found in ref. [158) . 
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In the continuum limit, describing the physics at length scales much larger than £ e , this random 
walk may be described as a Brownian motion. It can be shown that (Act) is related to the 
time integrated probability for an electron to come back to its initial position (see ref. jHSl f° r a 
heuristic derivation or ref. |I] for a recent derivation in space representation) : 

9 p 2 j~> poo 

(Acx) = — — dtV(f,t\r,0)e- t/T * . (2) 

where D is the diffusion constant and the factor 2 accounts for spin degeneracy. V(r, t\f ', 0) is 
the Green's function of the diffusion equation 

(jt-DAj V(f,t\r',0)=5(f-r')5(t) . (3) 

In eq. © the exponential damping at large time describes the lack of phase coherence due to 
inelastic processes. The phase coherence time can be related to the phase coherence length 
by = Dt v . We have also introduced a cut-off at short time in eq. ©, that takes into account 
the fact that the diffusion approximation is only valid for times larger than r e . 

2.2 Planar Brownian motion : stochastic area, winding and magnetoconduc- 
tance 

Weak localization corrections are directly related to the behaviour of the probability of return 
to the origin and thus to recurrence properties of the diffusion process. It therefore follows that 
dimension d = 2 plays a very special role. Consider for instance a thin film whose thickness 
a is much less than L v . The sample is effectively two dimensional and the correction to the 
conductivity is given by 

(Aa) = -^ h HL v /Q. (4) 

On probabilistic grounds the logarithmic scaling is of course not unexpected here, it is the 
same logarithm which occurs in asymptotic laws of the planar Brownian motion |181l I151j . The 
neighbourhood recurrence of the planar Brownian motion favors the quantum interference effects 
and leads to a reduction of the electrical conductivity. A more striking effect is predicted if one 
applies a constant and homogeneous magnetic field B over the sample. In this case the classical 
action contains a coupling to the magnetic field S = eBA where A is a functional of the path 
given by the line integral 

A = \ j ( xd y ~ y^ x ) ■ ( 5 ) 

Properly interpreted, this line integral is nothing but the stochastic area of planar Brownian 
motion, whose distribution was first computed by P. Levy before the discovery of the Feynman 
path integral. The weak localization correction reads 

(M«)) - No)) = I j e ~ t/Tv v 1 " E [^ eBA/h ] ) ■ (6) 

The coupling to the magnetic field now appears with an additional factor 2 coming from the fact 
that the weak localization describes quantum interferences of reversed paths. The expectation 
E[- • • ] is taken over Brownian loops for a time t. The magnetoconductivity is given by |12 [ l34 | l55] 



(Aa(B)) - <M0)) ' 



2 



2ir 2 h 



4e 2 (BLl 



3 h 



(7) 
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where ip{z) is the digamma fonction and 4>q = h/e is the quantum flux. The r.h.s corresponds to 
the weak field limit B <C (fto/L^. This increase of the conductivity with magnetic field is opposite 
to the behaviour expected classically. This phenomenon is called "positive (or anomalous) 
magnetoconductance" . Experimentally this effect is of primary importance : the magnetic field 
dependence allows one to distinguish the weak localization correction from other contributions 
and permits to extract the phase coherent length L v . This expression fits the experimental 
results remarkably well 7 |34j . 

Strictly speaking eq. only holds for an infinite sample. In the case of bounded domains 
the variance depends on the geometry of the system and can be computed explictly for rectangles 
and strips [20] • The case of inhomogeneous magnetic fields can be treated along the same line. 
Consider for instance a magnetic vortex carrying a flux <f> and threading the sample at point 
0. In this case the functional which is involved is not the stochastic area but the index of the 
Brownian loop with respect to (winding number). The probability distribution of the index 
has been computed independently by Edwards [HS] m the context of polymer physics and Yor 
|227| in relation with the Hartman- Watson distribution. For a planar Brownian motion started 
at r and conditioned to hit its starting point at time 1, the distribution of the index when n 
goes to infinity is 

Proba [Ind = n] ~ — ^ K (r 2 ) e" r2 , (8) 

where Ko(cc) is the modified Bessel function of second kind (MacDonald function). The fact 
that the even moments of this law are infinite is reflected in the non analytic behaviour of the 
conductivity 

<A(r(0)> - <A<t(0)> oc |0| . (9) 

It is interesting to compare with eq. (J3| where the quadratic behaviour is given by the sec- 
ond moment of the stochastic area 169 . The behaviour given by eq. © has been observed 
experimentally |107| 131 j and a theoretical interpretation is provided in ref. |182j . 



2.3 Dephasing due to electron-electron interaction and functionals of Brow- 
nian bridges 

Dephasing in a wire : relation with the area below a Brownian bridge.- Another interesting 
example of a non trivial connection between a physical quantity and a functional of the Brownian 
motion occurs in the study of electron-electron interaction and weak localization correction in 
a quasi-ld metallic wire. In ref. [Oj, APtshuler, Aronov & Khmel'nitzkh (AAK) proposed to 
model the effect of the interaction between an electron and its surrounding environment by the 
interaction with a fluctuating classical field. The starting point is a path integral representation 
of the probability V(r,t\f', 0) in which is included the effect of the fluctuating field that brings a 
random phase. After averaging over Gaussian fluctuations of the field, given by the fluctuation- 
dissipation theorem, AAK obtained 

/A \ n e2 -° Pj, -<vt r (t)=I r, I \ -377 Jn dritr) 2 -^/,,' dr \x(r)-x(t-r)\ M 

(Act) = -2—— / die 7 / Vx(r)e 4DJo y ' a o«w Jo 1 w v >\ ^ 

Jo Jx(0)=x 

where T is the temperature (the Planck and Boltzmann constants are set equal to unity h = 
ks = 1), and S w is the area of the section of the quasi-ld wire. In contrast with eq. © where 
the loss of phase coherence was described phenomenologically by an exponential damping 8 with 

7 Note that in materials with strong spin orbit scattering, like gold, the effect is reversed and a negative 
magnetoconductance is observed at small magnetic fields |127II34| . 

8 Note that the introduction of an exponential damping describes rigorously several effects : the loss of phase 
coherence due to spin-orbit scattering or spin- flip I127| . the penetration of a weak perpendicular magnetic field in 
a quasi Id wire In this latter case the effect of the magnetic field is taken into account through 7 = ^(eBW) 2 
where W is the width of the wire of rectangular section. 
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the parameter r„, here the electron-electron interaction affects the weak localization correction 
through the introduction in the action of the functional of the Brownian bridge 9 



At 



dr \x(t) — x(t 



(11) 



The additional damping exp — yt in eq. lfTU|) describes the loss of phase coherence due to other 
phase breaking mechanisms. By using a trick which makes the path integral local in time, AAK 
have computed explicitely the path integral and obtained : 



(Aa) 



L 



N 



Ai(7T Af ) 
Ai'(7rjv) 



(12) 



where Ai(z) is the Airy function. The Nyquist time tn = ( ^^ ) 2 ^ 3 , gives the time scale over 
which electron-electron interaction is effective and therefore plays the role of a phase coherence 
time. We have also introduced the corresponding length = \ / TJtn. Note that the T depen- 
dence of tn oc T -2 / 3 directly reflects the scaling of the area with time : At = i 3//2 -4.i. We stress 
that the AAK theory makes a quantitative prediction for the dependence of the phase coherence 
length L^v as a function of the temperature, which has been verified experimentally for a wide 
range of parameters : on metallic (Gold) wires jSlllTKnj (the behaviour of L N oc (S w /T) 1 / 3 was 
observed before in Aluminium and Silver wires |224j ). and on wires etched at the interface of 
two semiconductors |213| . 





Figure 1: Left : The functional At gives the absolute area below a Brownian bridge x(t). Right : 
The functional At measures the area between the Brownian bridge x(t) and its time reversed 
counterpart x(t — r). The two functionals are equal in law [eq. H4\ )I- 

It is interesting to point out that the result of AAK can be interpreted as the Laplace 
transform of the distribution of the functional E[exp — pAt]- This functional represents the 
area between a Brownian bridge and its time reversed counterpart (c/. figure E^), where E[- • • ] 
describes averaging over Brownian bridges. The conjugate parameter p is played in by the 
temperature : p = ■§^§~- The result ()12|) has also been derived in the probability literature. Let 
us consider the functional 

At= I dr\x(r)\ , (13) 
Jo 

giving the absolute area below the Brownian bridge starting from the origin, x(0) = x(t) = (c/. 
figure^). The double Laplace transform of the distribution of At has been computed first by 
Cifarelli and Regazzini [SH] in the context of economy and independently by Shepp 195,, and led 
to / °° ^e-^E[e-^ At ] = -v^Xj^T), wnich is equivalent 10 to the result (JUJ). The connection 

9 A Brownian bridge, (x(t), ^ r ^ 1 x(0) = x(l)), is a Brownian path conditioned to return to its starting 
point. 

10 Note that the in the integral computed by Shepp corrects the fact that the averaging over Brownian 

curves in HOH is not normalized to unity, whereas E[- ■ ■ ] is. 
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between the two results is clear from the equality in law 



(law) 



A t 



(14) 



which follows from 



x(t) — x(t — t) 



(law) 



x(2t) for rG[0,(/2], 



(15) 



The relation (j!4j) allows us to understand more deeply the trick used by AAK to make the path 
integral (|10[) local in time. The distribution of the area At has been also studied by Rice in 
ref. ^83 . Interestingly, the inverse Laplace transform of eq. (|12|) obtained by Rice has been 
rederived recently independently in ref. |168j in order to analyze the loss of phase coherence due 
to electron-electron interaction in a time representation. 

The study of statistical properties of the absolute area below a Brownian motion was first 
addressed by Kac |135j for a free Brownian motion, long before the case of the Brownian bridge 
studied by Cifarelli & Regazzini and by Shepp. Later on, it has been extended to other func- 
tionals of excursions and meanders 13 |133[ I179j , which arise in a number of seemingly unrelated 
problems such as computer science, graph theory [HSj and statistical physics 1601 116 J] . 

Dephasing in a ring. - Very recently, the question of dephasing due to electron-electron interac- 
tion in a weakly disordered metal has been raised again in ref. |159j . In particular it has been 
shown that the effect of the geometry of the ring and the effect of electron-electron interaction 
combine in a nontrivial way leading to a behaviour of the harmonics of the magnetoconductance 
that differs from the one predicted by eq. (J2J in ref. • For the ring, the functional describing 
the effect of electron-electron interaction on weak localization is now given by 



Jo \ ^ / 

instead of the area defined by eq. Q13JI . x(t) is a Brownian path 14 on a ring of perimeter L such 
that x(0) = x(t) = 0. The question has been re-examined in more detail in ref. |211| . where the 

double Laplace transform of the distribution J °° dt e~ 7 * "^e - { » E^e -7 ^] has been derived. 
E n [- • •] denotes averaging over Brownian bridges defined on a circle, with winding number n. 
Note that the Laplace transform of the distribution E n [e -7 ^*] has been also studied in ref. |211j . 

3 Exponential functionals of Brownian motion and Wigner time 
delay 

3.1 Historical perspective 

Exponential functionals of the form, 



11 For a given process x(t) the two functionals At and At are obviously different, however they are distributed 
according to the same probability distribution. They are said "equal in law" . 

12 The relation 11151 is easily proved by using that a Brownian bridge (x(r), ^ r ^ 1 1 £c(0) = £c(l) = 0) can be 
written in terms of a free Brownian motion (B(t), t ^ ] B(0) — 0) as : x(t) = B(t) — r B(l). 

13 An excursion is a part of a Brownian path between two consecutive zeros and a meander is the part of the 
path after the last zero. 

14 As for the case of the wire, the functional involved in the study of dephasing in the ring |211| involves the 
difference x(t) — x(t — t) instead of the bridge x(t). We have used the equality in law 115H which implies the 
following property : given a Brownian bridge (x(r), ^ r ^ t \ x(0) = x(t) = 0), for any even function f(x) we 




(16) 




(17) 



have J*dTf(x(r)-x(t-T)) = / 'dr/(a;(r)). 
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where (B(s), s ^ 0, -B(O) = 0) is an ordinary Brownian motion have been the object of many 
studies in mathematics [228| . mathematical finance and physics [HI]. in the- physics of classical 
disordered systems, the starting point was the analysis of the series (Kesten variable) 

Z = z\ + z\z 2 + ziz 2 z 3 H (18) 

where the Z\ are independent and identical random variables. Several papers have been devoted 
to the study of this random variable in the physics |68M71| and mathematics 113811217] literature. 
It was realized that Z is the discretized version of A { £ defined in eq. (|T7j) . which may be 
interpreted as a trapping time in the context of classical diffusion in a random medium. The 
tail of the probability distribution P(Z) controls the anomalous diffusive behaviour of a particle 
moving in a one-dimensional random force field [15] . The functional also arises in the study 
of the transport properties of disordered samples of finite length |170[ 1173] . In the context of 
one dimensional localization the fact that the norm of the wave function is distributed 15 as 
is mentioned in the book of Lifshitz et al |157j . 



3.2 Wigner time delay 

As discussed in the introduction, the localization of quantum states in one dimension is well 
understood. However since real systems are not infinite, asking about the nature of the states 
which are not in the bulk is a perfectly legitimate question. It was pointed out by Azbel in 
ref. [U that, when a disordered region is connected to a region free of disorder, the localized 
states acquire a finite lifetime which shows up in transport through sharp resonances, refered as 
"Azbel resonances" . This picture was used in ref. [T32] where it was argued that these resonances 
could be probed in scattering experiments and lead to an energy dependent random time delay 
of the incident electronic wave. Motivated by this result and by developments in random matrix 
theory we consider the scattering problem for the one dimensional Schrodinger equation (in units 
h = 2m = 1) 

d 2 

— r^Mx) + V(x)ip(x) = Eibix) (19) 
dx z 

defined on the half line x ^ with the Dirichlet boundary conditions V'(O) = 0. The potential 
V(x) has its support on the interval [0, L]. Outside of this interval, the scattering state of energy 
E = k 2 is given by 

iPe(x) = -L= L-^ x - L ) + e iH*-L)+is\ ^ (20) 
yhvE ^ ' 

where h = 2tt is the Planck constant (in unit h = 1) and ve = dE/dk = 2k the group velocity. 
Eq. (|2U|) represents the superposition of an incoming plane wave incident from the right and a 
reflected plane wave caracterized by its phase shift 5(E). The Wigner time delay, defined by the 
relation r = d5/dE, can be understood as the time spent by the wave packet of energy E in the 
disordered region (this interpretation is only valid at high energy ; see refs. |125| I50| I149| 169] 
for review articles on time delay and traversal times). This representation of the time delay has 
been used in many papers both in the mathematics and physics literature. Starting from this 
representation one can derive a system of stochastic differential equations which can be studied 
in certain limiting cases. However, to understand the universality of the statistical properties of 
the Wigner time in the high energy limit, we follow a different approach below. 

15 When fi > 0, possesses a limit distribution for t — > oo. Moreover, we have the equality in law (see 

footnote[TT]| : A ( £ } = 1/7 (ai) , where 7 (m) is distributed according to a T-law : P(j) = r7^7 M_1 e" 7 . 
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Universality of time delay distribution at high energy 



Using the expression of the scattering state (|20j) and the definition r = d5/dE, we can obtain 
the so-called Smith formula |98l I196j relating the time delay to the wave function in the bulk : 

r(E) = 2n J dx|^(x)| 2 -± S m8(E) . (21) 

Following ref. ^3] one can parametrize the wave function in the bulk in terms of its phase and 
modulus. For this purpose we rewrite the Schrodinger equation as a set of two coupled first order 
differential equations, dipE/dx = ip' E and dip' E /dx = (V(x) — E)ipE, and perform the change of 
variables : iPe(x) = e^ x ' sm8(x) and iP'e( x ) = ke^ x ' cos 6{x). One obtains a new set of first 
order differential equations d6/dx = k — ^V(x)s'm 2 9 and d^/dx = j^V(x) sin29. Up till now 
everything is exact. If we now consider the high energy limit 16 we can neglect the second term 
on the r.h.s of eq. ()21j) and integrate out the phase which is a fast variable ; we obtain 

r = - f L dxe 2 ^-^. (22) 
k Jo 

This representation of the time delay holds for any realization of the disordered potential. More- 
over, one can prove under rather mild conditions on the correlations of the random potential 
that £(x) is a Brownian motion with drift £(x) = x/X + a/1/A B(x) where B(x) is an ordinary 
Brownian motion and A the localization length. Using the scaling properties of the Brownian 
motion gives the following identity in law 



r 



(law) 



H iA d,e-^)«> = ^™ , (23) 



where has been defined in eq. (|17|). This representation of the time delay as an exponential 
functional of the Brownian motion, first established in ref. by a different method, allows us 
to obtain a number of interesting results |65l 207 : 

(i) Existence of a limit distribution for fixed 17 r and L — > oo 

P(r) = ^A_e-^. (24) 

This result is reminiscent of the random matrix theory prediction in spite of the fact that this 
theory does not apply to systems that are strictly one-dimensional 18 . 

(m) Linear divergence of the first moment, (r) = L/k, and exponential divergence of the higher 
moments (r n ) oc Q 2n i n ~ l ) L /^ . This divergence reflects a log-normal tail of the distribution for a 
finite length L. 



Time delay and density of states 

The relation (J21[) also provides another interpretation to these results. In the situation under 
consideration, the scattering state is directly related to the local density of states (LDoS) by 
p(x;E) = (x\5(E — H)\x) = |^(x)| 2 . Therefore the time delay can be interpreted as the DoS 

16 To define precisely the high energy limit, let us introduce the integral of the correlation function of the disorder 
w = J dx (V(x)V(O)). The high energy limit corresponds to k ^> w 1 ^ 3 . 

17 From the remark of the footnote 1151 we notice that 1/r is distributed according to an exponential law. 

18 The random matrix theory describes the regime of weak localization or systems whose classical dynamics 
is chaotic. The distribution of time delay in this framework has the same functional form P(r) oc Q~ T °^ T 
[1021 11141 1171 1174| with two differences : (i) the exponent 2 + fi > 2, (it) the time scale to- 
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of the disordered region 19 . This establishes a relation between the results given above and the 
work of Al'tshuler & Prigodin ^1] where the distribution of the LDoS was studied for a white 
noise potential using the method of Berezinskh 20 . 

Time delay for Dirac Hamiltonian at the threshold energy 

The representation (|23|) of the time delay as an exponential functional of the Brownian motion 
only holds in the weak disorder (i.e. high energy) limit. A similar representation can be derived 
for the random mass Dirac model at the middle of the spectrum. The Dirac Hamiltonian is 

H D = a 2 i^+a 1 ^(x) , (25) 

where a-i are the Pauli matrices. (p(x) can be interpreted as a mass 21 . The Dirac equation 
H-Qip = kip possesses a particle-hole symmetry reflected in the symmetry of the spectrum with 
respect to k = (the dispersion relation of the Dirac equation is linear in absence of the mass 
term and energy is equal to momentum). The middle of the spectrum is an interesting point 
where the divergence in the localization length signals the existence of a delocalized state (see 
footnote I28j) . These properties should therefore show up in the probability distribution of the 
Wigner time delay. Since the dispersion relation of the Dirac equation is linear, the time delay 
is now defined as r = d<5/d/c. The following representation has been derived in ref. ^ 



T 



L 



2 / dxe 2 ^^ y)dy . (26) 



o 



In contrast with eqs. (|22l2^j) . which have been obtained after averaging over the fast phase 
variable, the representation l)26|) is exact (for k = 0). If (j)(x) is a white noise, then 

In this case since the drift vanishes, it is known |17U1 1228] that there is no limiting distribution 
when L — > oo. The physical meaning of the lognormal tail 

in the limit r — > oo is discussed in detail in ref. 198 . 



4 Extreme value spectral statistics 

The density of states (DoS) and the localization length of a one-dimensional Hamiltonian with a 
random potential can be studied by the phase formalism |151 1157) (see also ref. |158j for discrete 

19 The relation between time delay and DoS is sometimes refered as Krein-Friedel relation |97l I144| . This 
relation has been originally introduced in ref. |37| in the context of statistical physics (see also ref. |67| and §77 
of ref. 148 ). Note that it has been recently discussed in the context of graphs in refs. 205 206 208 . 

20 The method of Berezinskii blocks has been introduced to study specifically the case of white noise disordered 
potential in one dimension ;32^ . This powerful approach has been widely used and has allowed to derive numbers 
of important results like in refs. [1121 111.5] for example. 

21 The Hamiltonian 125H with 4>(x) a white noise was introduced by Ovchinnikov & Erikmann |176| as a model 
of one- dimensional semiconductor with a narrow fluctuating gap. It is interesting to point out that the Dirac 
equation Huip = kip also has an interpretation in the context of superconductivity, as linearized Bogoliubov- 
de Gennes equations for a real random superconducting gap 4>(x). Finally it is worth mentioning that more 
general Id Dirac Hamiltonians with several kinds of disorder (mass term, potential term and magnetic field) have 
been studied in refs. |38l 139) . 
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models). In the previous section we pointed out that the physics of the strong localization in 
one-dimension can also be probed from a different angle by considering the scattering of a plane 
wave by the random potential. This has led to the study of the Wigner time delay. In this 
section we show that the phase formalism also allows describing finer properties of the spectrum 
such as the probability distribution of the n-th eigenvalue. This is an instance of an extreme 
value statistics : given a ranked sequence of ftf random variables x\ ^ 22 ^ • • • ^ %Af, the 
problem is to find the distribution of the n-th of these variables in a given interval. In the 
particular case of uncorrelated random variables, the extreme value distributions have been 
studied by E. Gumbel |119M120"lll21| . This problem becomes much more complicated when the 
random variables are correlated, a case which has recently attracted a lot of attention. Such 
extreme value problems have appeared recently in a variety of problems ranging from disordered 
systems |44[ I54[ I70| to certain computer science problems such as growing search trees |162| . 
In the case of the spectrum of a random Hamiltonian, the eigenvalues are in general correlated 
variables, apart from the case of strongly localized eigenstates |167j . Below we show that, in 
the one dimensional case, this problem is related to studying the first exit time distribution of 
a one-dimensional diffusion process. 



4.1 Distribution of the n-th eigenvalue : relation with a first exit time dis- 
tribution 

We consider a Schrodinger equation Hip(x) = E(p(x) on a finite interval [0, L]. The spectral 
(Sturm-Liouville) problem is further defined by imposing suitable boundary conditions. We 
choose the Dirichlet boundary conditions c/?(0) = f(L) = 0. The spectrum of H is denoted by 
Spec(H) = {Eq < Ei < Ei < • • • }. Our purpose is to compute the probability 

W n (E) = (5(E - E n )} (29) 

for the eigenvalue E n to be at energy E (the bracket (• • •) means averaging over the random 
potential). Note that the sum of these distributions j^^ n W n (E) = p(E) is the average DoS 
per unit length. We now show how the calculation of W n {E) can be cast into a first exit time 
problem. 



Random Schrodinger Hamiltonian 

We consider the Hamiltonian 

H = --f^ + V(x), (30) 
ax A 

where V{x) is a Gaussian white noise random potential : (V(x)) = and (V(x)V(x')) = 
w S(x — x'). 

We replace the Sturm-Liouville problem by a Cauchy problem : let ip(x;E) be the solution 
of the Schrodinger equation Htf)(x;E) = Eip(x;E) with the boundary conditions vp(0;E) = 
and 4^ip(0;E) = 1. The boundary condition ip(L;E) = is fulfilled whenever the energy E 
coincides with an eigenvalue E n of the Hamiltonian. In this case, the wave function (p n (x) = 

r L i 1/2 

tp(x; E n )/ f Q dx' tp(x'; E n ) 2 has n nodes in the interval ]0, L[, and two nodes at the bound- 
aries. 

Let us denote by £ m (m ^ 1), the length between two consecutive nodes. We consider the 
Ricatti variable 

z(x;E) = ^-ln\^x;E)\, (31) 
which obeys the following equation: 

^-z = -E - z 2 + V(x) , (32) 
dx 
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Figure 2: The probability for the n-th level E n of the spectrum to be at E is also the probability 
for the n + 1-th node of ' ip(x; E) to be at L. 

with initial condition z(0;E) = +oo. This equation may be viewed as a Langevin equation for 
a particle located at z submitted to a force —dU(z)/dz deriving from the unbounded potential 

U(z) = Ez + j (33) 

and to a random white noise V{x). 

Each node of the wave function corresponds to \z{x)\ = oo. At "time" x = the "particle" 
starts from z(0) = +oo and eventually ends at z{l\ — + ) = — oo after a "time" t\. Just after 
the first node it then starts again from z{l\ + + ) = +oo, due to the continuity of the wave 
function. It follows from this picture that the distance l m between two consecutive nodes may 
be viewed as the "time" needed by the particle to go through the interval ] — oo,+oo[ (the 
"particle" is emitted from z = +oo at initial "time" and absorbed when it reaches z = —oo). 
The distances £ m are random variables, interpreted as times needed by the process z to go 
from +oo to — oo. These random variables are statistically independent because each time the 
variable z reaches — oo, it loses the memory of its earlier history since it is brought back to the 
same initial condition and V(x) is <5-correlated. This remark is a crucial point for the derivation 
of W n (E) . Interest in these random "times" i m lies in their relation with the distribution of the 
eigenvalues. Indeed the probability that the energy E n of the n-th excited state is at E is also 
the probability that the sum of the n + 1 distances between the nodes is equal to the length 
of the system: L = Yl^ri^i^m (this is illustrated on figure HJ). Since the t m are independent 
and identically distributed random variables Prob L = Y^=i is readily obtained from the 

distribution P(i) of one of these variables. 

We introduce the intermediate variable C(z) giving the time needed by the process starting 
at z to reach — oo : therefore we have I = £(+oo). The Laplace transform of the distribution 
h{a,z) = (e~ aC | z{0) = z; z(C) = -oo) obeys [Ml : 

G z h(a, z) = a h(a, z) , (34) 

where the backward Fokker-Planck generator is 

G z = -U'(z)d z + |a 2 2 . (35) 
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The boundary conditions are d z h(a, z)\ z=+OQ = and h(a, —oo) = 1. By Laplace inversion we 
can derive P{£) since /i(a,+oo) = d£ P(£)e~ ai . 

The solution of this problem is given in ref. |2U4| . Here we only consider the limit E — > — oo, 
which corresponds to the bottom of the spectrum. In this regime the dynamics of the Ricatti 
variable z can be read off from the shape of the potential (see figure EJ). 





\U(z) . 












E<0 



Figure 3: The potential 4.V.V)) related to the deterministic force to which the Ricatti variable is 
submitted in eq. \3ty) . 

The particle falls rapidly in the well from which it can only escape by a fluctuation of the 
random force. Therefore the time needed to go from +oo to — oo is dominated by the time spent 
in the well, which is distributed according to the Arrhenius formula. It follows that the average 
"time" , which is related to the inverse of the integrated density of states (IDoS) per unit length 
N(E) = J^dE' p(E'), behaves as {£) = N(E)- 1 ~ exp ^(-£) 3 / 2 . This simple picture, 
first provided by Jona-Lasinio |134| allows recovering the exponential tail obtained by several 
methods in refs. |991 11241 1151 11571 H3Uj 22 . The distribution of the length £ is a Poisson law (see 
jl()5j or appendix of ref. |2()4| ): 

P{£) = N{E) exp -£ N(E) . (36) 
Using this result we can show that ^204 

W n (E) = Lp(E) (LiY( f }) " e~ LN ^ . (37) 
re! 

This result has a clear meaning : Lp(E) gives the probability to find any level at E and the 
factor ^-e _x "compels" the number of states below E, x = LN(E), to be close to n. We may 
go further and write 

W ^ E) = JK (-k*-) ' (38) 

where the typical value of the energy is E^ P (L) = -(^lnL) 2 / 3 , while the scale of the fluctua- 
tions reads SE n = w . = (3 In L) -1 / 3 , where L = ^r^TTv ^ ne function 



(n + 1) 



uj n (X) = v ^\> * exp( v / nTTX-(n + l)e x /^+ T ) (39) 



has the form of a Gumbel law for uncorrelated random variables. The fact that the eigenvalues 
are uncorrelated is a consequence of the strong localization of the eigenfunctions |167j . 

22 It is interesting to quote ref. [145] where it was shown that this non perturbative exponential tail can be 
obtained by a simple counting method of "skeleton" diagrams assuming that they all have the same value (this 
approach was also applied to the 3d case). 
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The work summarized in the above paragraph, which appeared in ref. |2(J4| . generalizes the 
result of McKean |165j for the ground state. 

A similar result was obtained by Grenkova et al |117j for the model of <5-impurities with 
random positions, in the limit of low impurity density which has no counterpart in the model 
we considered here 23 . However both models describe the same physics of strongly localized 
eigenstates. 

Supersymmetric random Hamiltonian 

Random Schrodinger operators have been investigated through a wide range of models. De- 
pending on the physical context, there are indeed many ways to model the disorder. In the 
previous paragraph we have assumed that the potential V{x) is a white noise. We have also 
mentioned that if the potential is a superposition of (5-potentials randomly distributed along the 
line, although the spectra of the two models are quite different 24 , the extreme spectral statistics 
are the same. The localization properties and the statistics of the time delay at high energy 
are also similar for both models |203| I207| . These two models belong to the same class of ran- 
dom Hamiltonians with a random scalar potential with short range correlations. They are both 
continuous versions of discrete tight binding models with on-site random potential. This is the 
case of so-called diagonal disorder, since the random potential appears on the diagonal matrix 
elements of the tight binding Hamiltonian in the basis of localized orbitals 25 . 

Other interesting models can be constructed by introducing disorder in the hoppings, instead. 
We refer to such models as off-diagonal disorder. The Dirac Hamiltonian (|25|) introduced above 
provides a continuum limit of such a model 26 . Dirac Hamiltonians appear naturally in several 
contexts of condensed matter physics. The existence of symmetries in the Dirac Hamiltonian 27 
(particle-hole, chiral,...) can lead to interesting features in the presence of disorder, which 
has attracted some attention (see for example refs. |1721 |7J 0S1 H33 E3 HUE EB]> the last of 
these references gives a brief overview). The square of the Dirac Hamiltonian, H 2 ^ = — d 2 + 
4>{x) 2 + o~s(f)'(x), is related to the pair of supersymmetric Schrodinger isospectral Hamiltonians 
H± = — d 2 + (j)(x) 2 ± (j)'(x). When (4>(x)) = we may forget the sign and simply consider : 



Such Hamiltonians appear in a variety of problems ranging from the Id classical diffusion in a 
random force |421ll28j . electronic structure of polyacetylene |2UU| . spin Peierls chains |91U92lH99| 

23 See footnote El 

24 The model of <5-impurities with random positions was introduced and studied by Schmidt 188 but often 
refered as the Frisch & Lloyd model | 99| . Compared to the model for white noise potential, which is characterized 
by one parameter, the Frisch & Lloyd model is characterized by two parameters : the strength of the 5-potential, 
and their density. In the limit of high density of impurities this model is equivalent to the white noise potential 
model. The limit of low density presents different spectral singularities (Lifshitz singularity and, for negative 
weight of 5-potentials, an additional Halperin singularity in the negative part of the spectrum). 

25 Note that the continuum limit of tight binding Hamiltonian with diagonal disorder can lead to different 
continous models. Let us consider the discrete model Hi,j = —Sij+i — Si,j-i + Si.jVi. For Vi = the spectrum is 
Ek = — 2cos(fc), with fc G] — 7r, n]. (A) If the continuous limit is taken by considering the band edge (k ~ 0), one 
is led to the continuous model 13011 . (B) If the band center is considered instead (k ~ 7r/2), the spectrum can be 
linearized and one is led to a Dirac Hamiltonian, like in refs. |137U17r)] . This point has been recently rediscussed 
in ref. (155) . 

26 The continuum limit of a tight binding Hamiltonian with random hoppings has been discussed in ref. |2(JU| 
(see also the review in ref. |169| L As pointed in footnote 1251 random Dirac Hamiltonians can also appear as 
continuum limit of the band center of a discrete models with diagonal disorder. 

27 A complete classification of symmetries of disordered Hamiltonians extending the famous Wigner-Dyson 
ensembles of random matrix theory has been provided in refs. |216l 12291 |7j (see also the recent review article 




(40) 
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and also in a continuum limit of the random field XY model 122 (see ref. |66j for a short review 
on supersymmetric disordered quantum mechanics). Spectral and localization properties have 
been studied in detail when (f>(x) is a white noise [17(51 11571 1411 142j and also when (f>(x) has a 
finite correlation length as in a random telegraph process jSJ (this last case has found some 
application in the context of spin Peierls chains). 

In the high energy limit, the localization properties and the statistics of the time delay 
do not show any difference with respect to the case of diagonal disorder. However the low 
energy properties are quite different. This is easily understood by noticing that, due to its 
relation to the Dirac Hamiltonian, the supersymmetric Hamiltonian can be factorized as Hg = 
Q^Q where Q = —d x + <p{x) and = d x + 4>(x). In particular, such a structure enforces a 
positive spectrum : Spec(-ffs) C M + . When <j){x) is white noise of zero mean, (<j)(x)) = and 
(4>(x)(J)(x')) = gd(x — x'), the case considered below, the DoS presents a logarithmic Dyson 
singularity : the integrated density of states (IDoS) reads N(E) ~ 1/ln 2 E |1761 l4*Tl E2] . while 
the localization length diverges logarithmically 28 X(E) ~ ln(l/E). 

It is interesting to investigate the extreme value statistics of the spectrum in the low energy 
regime, where we expect properties quite different from the one obtained for the diagonal dis- 
order. The derivation follows closely that in the previous case, however the relevant random 
processes and the approximations are different. The first step is to decouple the Schrodinger 
equation Hsf{x) = k 2 (p(x) into two first order differential equations {i.e. go back to the Dirac 
equation) : 

Qt x ( x ) = k<p(x) (41) 
Q<p(x) = k x (x). (42) 

Then we may use the phase formalism by introducing a phase variable and an envelope variable : 
<p(x) = e^ x ' sm , &(x) and xi x ) = ~ cos , &(x). The phase variable obeys a stochastic differ- 
ential equation with a noise multiplying a trigonometric function of the phase. For convenience 
we introduce an additive process ((x) defined as £(a?) = ±^ln| tan$(x)|. The sign + is chosen 
for ($ mod n) G [0,7r/2] and the sign — for (i? mod tt) £ [tt/2,tt]. This new process obeys the 
stochastic differential equation 

— C = A; cosh 2C ± 0(s) • (43) 
dx 

Between two nodes of the wave function, the variable £ crosses twice the interval ] — oo,+oo[. 
Note that when <fi(x) is a white noise of zero mean, the sign ± can be disregarded. The study 
of the time required by the process to cross the interval can be performed by the same method 
as above : we introduce the "time" A(£) needed to go from £ to +oo. The Laplace trans- 
form of the distribution h(a,() = (e~ aA | C(0) = £; C(A) = +oo) obeys a diffusion equation 
(fccosh2£ + |<9 2 )/i(a,£) = a h(a,() that involves the backward Fokker-Planck operator re- 
lated to the Langevin equation (|43|). 

In the low energy limit k <C g we expect that most of the "time" A = A(— oo) is spent in 
the region where the potential is almost flat (see figure 0J) , therefore we replace the diffusion 
equation for h(a,Q by the free diffusion equation |9 2 /i(a,C) = ot /i(a,C) on the finite interval 
( £ [C-,C+]; with a reflecting boundary condition at one side d^h(a, £-) = and h(a, C+) = 1 
at the other side (which corresponds actually to the absorption at £ = £ + ). The coordinates C± 
are the points where the deterministic force and the white noise have equal strengths. 

28 A more complete picture of the properties of this model at E — is given in several works : (A) moments 
and correlations of the zero mode are studied in refs. [1941 (B) Additionally to the statistical properties of 
the time delay , eq. (1271 . (C) the distribution of the transmission probability is derived in ref. |198| . In particular 
it was shown that the average transmission through an interval of length L decreases like (T) oc 1/ \fh, which is 
slower than in the diffusive regime. (D) The existence of a finite conductivity was demonstrated in ref. 
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Now, we can obtain h(a, £) straightforwardly. 



cosh y/a/N(E) ' 



(44) 



where N(E) = g/2ln 2 (g/k) is the IDoS per unit length. 



mo 



Hg/k) 
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Figure 4: The potential related to the deterministic force in eq. 14$) felt by the variable £. 

If A = A(— oo) is the "time" needed to cross the interval, the distance between two nodes of 
the wave function is a sum of two such (independent) random variables £ = A\ + A2. Therefore 
its distribution is given by inverse Laplace transformation of (e~ ai ) = 1/ cosh 2 y / a/N(E). We 
obtain : 

P(i) = N(E)w (N(E)£), (45) 
where wq{x) is the inverse Laplace transform of cosh -2 yfs : 



Y(x) J2 [vr 2 (2m+ 1) 



m=0 



4 Y(x) 



£(-D 



m+1 m 2 e _m ' x 



2] Q-Wm+lfx 

4 Y(x 



~ vr xe 4 



m=l 



x^O \[k ^ 2 



where Y(x) is the Heaviside function. 

We find the distribution of the ground state energy : 



W {E) = L P (E)vdq(LN(E)) 
with the following limiting behaviour : 



W (E) ~ 



8vr 2 <7 2 L 2 



1 in 2 (g 2 /E) 
exp - 



Eln 5 {g 2 /E) GXP 21n 2 (<? 2 /£) 



2gL 



for £ <C g 2 e~^ L 

for g 2 e-^ Z < £ < # 2 



(46) 
(47) 



(48) 

(49) 
(50) 



In ref. |2U4j . an integral representation of I4 / n (£') is also given together with an explicit expression 
for W X (E). 

It is also interesting to point out that the ground state distribution is characterized by 
a typical value Eq p ~ g 2 e~ gL , a median value EQ lcd ~ g 2 e~ v ^, and a mean value (Eq) ~ 
g 2 (gL) 1 / 2 e~ c (9 L ) 1 ^ 3 where C is a numerical constant (see ref. [US]). This latter expression has 
also been obtained in ref. |171j where upper and lower bounds were found using a perturbative 
expression for the ground state energy as a functional of 4>(x). 
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The distribution was obtained in ref. |15L)j in the context of the classical diffusion by 
using a real space renormalization group method. In this case the distribution is interpreted as 
the distribution of the smallest relaxation time. 

In summary, these two examples illustrate the fact that extreme value spectral statistics pro- 
vides an information on correlations of eigenvalues : the extreme value statistics of independent 
and identically distributed random variables have been classified by Gumbel |12L)j . Therefore ex- 
treme value distribution for eigenvalues that differ from one of the three Gumbel's laws indicates 
level correlations. 

5 Trace formulae, spectral determinant and diffusion on graphs 

5.1 Introduction 

Up to now we have discussed several questions related to the physics of weak and strong local- 
ization. In the previous section we have considered spectral properties while in the sections 
and |31 we mostly discussed transport properties. 

In the present section we are going to review several results related to the study of the 
Laplace operator on metric graphs. An object at the core of our discussion is the spectral 
determinant of the Laplace operator, formally defined as S(j) = det(7— A), where 7 is a spectral 
parameter. This quantity encodes the information on the spectrum of the Laplace operator on 
the graph. The subsection 15.21 recalls the basic conventions required to describe metric graphs. 
The subsections 15.31 and 15.41 review general results on trace formulae and spectral determinants. 
These subsections will appear at first sight quite technical and unrelated to the previous sections. 
However we will see that there exists a close relation with the question of quantum transport : 
the Laplace operator is the generator of the diffusion on the graph and its spectral properties 
play a central role in the study of transport. This connection, already evoked in the introduction 
and section [21 will be emphasized again in subsection 15.61 where the explicit relation between 
quantum transport and spectral determinant is recalled. 

One interest for spectral determinants is that they can be used as generating functions for 
various quantities characterizing the diffusion on the graph. This will be illustrated by using the 
connection with trace formulae (section l5.3j) and further exploited in subsection l5.6l in the context 
of quantum transport. The efficiency of the method stems from the fact that, even though £(7) 
seems to be a complicated object at first sight, involving an infinite number of eigenvalues, it can 
be expressed as the determinant of a finite size matrix. This relation, established by Pascaud 

6 Montambaux |177|I178] . allows computing easily and systematically the spectral determinant 
for arbitrary graphs. 

The study of the Laplace operator on metric graphs (or quantum graphs) is not restricted 
to transport and appears in many physical contexts ranging from organic molecules |187j . su- 
perconducting networks [B], phase coherent transport in networks of weakly disordered wires 
[ED HZH1 01 EH3, transport in mesoscopic networks [Ml [Ml E21 DH3 EDI EES1 M or metal- 
lic agregates [53| . We are not going to review this history and refer the interested reader to 
the references [HI I59( I141( 01 1143j (for a recent review see the appendix of Pavel Exner in 
ref. 0). Mathematical aspects of quantum graphs are discussed in the recent issues of Waves 
Random Media 14 (2004) and Journal of Physics A 38 (n°22) (June 2005). See in particular 
the refs. [TIHllTTTj . 

5.2 Description of metric graphs 

We collect some background material from the theory of graphs and establish our notations and 
conventions (illustrated on figure EJ. 
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Vertices, bonds, arcs.— Let us consider a network of B wires (bonds) connected at V vertices. 
The latter are labelled with Greek indices a, /?,... Therefore, the bonds of the graph can be 
denoted by a couple of Greek indices (a/3). The oriented bonds, denoted as arcs, will also play 
an important role. The two arcs related to the bond {a(5) will be labelled by a/3 and Pa or more 
simply with Roman letters i, j,... 

Adjacency matrix.— The basic object is the adjacency V x V matrix a a p characterizing the 
topology of the network : a a p = 1 if a and (3 are connected by a bond ; a a p = otherwise 29 . 
The connectivity of the vertex a, denoted by m a , is related to the adjacency matrix by m a = 
Yl/3 a af3- Summed over the remaining index, the adjacency matrix gives the number of arcs : 

Orbits.— A path is an ordered set of arcs such that the end of an arc coincides with the begining 
of the following arc. The equivalence class of all closed paths equivalent by cyclic permutations 
is called an orbit. An orbit is said to be primitive when it can not be decomposed as a repetition 
of a shorter orbit. 




Figure 5: An example of graph with 11 vertices and 13 bonds. The arrows show the orientations 
of the arc i and the reversed arc i. The blue curve is an example of (primitive) orbit with one 
backtracking. 

Scalar functions.— The graphs we consider here are not simply topological objects but have 
also some metric properties : the bond (af3) is characterized by its length l a p and identified 
with the interval [0, l a p] of R. A scalar function vb(x) is defined by its components ip a p{x a p) on 
each bond, where x a p is the coordinate measuring the distance along the bond from the vertex 
a (note that x a p + xp a = lap). 

Continuous boundary conditions.— When studying the Laplace operator A acting on scalar 
functions, boundary conditions at the vertices must be specified in order to ensure self adjointness 
of the operator. Let us introduce the notations : 

Paf3 = Paf3(x a/ 3 = 0) (51) 

^ = ^£(^ = 0), (52) 

dx a /3 

for the value of the function and its derivative at the vertex a, along the arc aj3. The continuous 
boundary conditions assume : 

(i) continuity of the function at each vertex : all the components (p a p for all (3 neighbours of a 
are equal. The value of the function at the vertex is denoted (f(a). 

(**) S/3 a ap t f' a p = A Q 9?(a) where the adjacency matrix in the sum constrains it to run over the 

29 Note that this definition assumes that two vertices are connected by at most one bond and that a bond never 
forms a loop. We emphasize that this does not imply any particular restriction on the topology of the graph : 
one can always introduce a vertex on a bond, which separates the bond into two bonds, without modifying the 
properties of the graph. The numbers of vertices and bonds can always be made arbitrary large and this is partly 
a matter of choice. The case of closed bonds (forming a loop) or vertices connected by multiple bonds requires a 
simple generalization of the formalism presented here (see for example ref. |177j and appendix C of ref. J2). This 
generalization allows in particular to minimize B and V, which makes the computation sometimes easier. 
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neighbouring vertices of a. Therefore the sum runs over all wires issuing from the vertex. The 
real parameter X a allows describing several boundary conditions : X a = oo enforces the function 
to vanish, <p(ct) = 0, and corresponds to the Dirichlet boundary condition. X a = corresponds 
to the Neumann boundary condition. We refer to the case of a finite X a as "mixed boundary 
conditions" . In the problem of classical diffusion, the Dirichlet condition describes the connection 
to a reservoir that absorbs particles, while the Neumann condition ensures conservation of the 
probability current and describes an internal vertex. 

General boundary conditions.— Thanks to the continuity hypothesis, the simple boundary 
conditions we have just described allow us to introduce vertex variables (f(ct) (this is convenient 
since the number of vertices V is usually smaller than the number of arcs 2B). However these are 
not the most general boundary conditions, and in general the different components associated 
to the arcs issuing from a given vertex do not have a priori the same limit at the vertex. The 
most general boundary conditions can be written as : 

C tp + D (p f = (53) 

where ip and ip' are the column vectors of dimension 2B collecting the variables (|51|) and (|52|). 
respectively. C and D are two square 2B x 2B matrices. In the arc formulation, the information 
on the topology of the graph is encoded in these two matrices. The self adjointness of the 
Laplace operator is ensured if these two matrices satisfy the following conditions |L-i9j : (i) CD^ 
is self-adjoint, (u) The 2B x AB matrix (C, D) must have maximal rank. Note that the choice 
of C and D is not unique. The introduction of boundary conditions without continuity at the 
vertices has been motivated physically in refs. [T321 HUHl 1521 ITS1 \M 1971] . 

Scattering theory interpretation.— We can give a more clear physical meaning to these 
conditions in a scattering setting. Let us consider the equation — A(p(x) = E(p(x) for a positive 
energy E = k 2 . We decompose the component on the bond as the superposition of an incoming 
and an outgoing plane wave : tpapixap) = I a f3^~ lkXal3 + O a @ e lkx <*p . It is convenient to collect 
the amplitudes in column vectors I and O. The vector / contains the incoming plane wave 
amplitudes and the vector O the outgoing plane wave amplitudes. Both are related by a vertex 
scattering matrix : O = Q I. The self- adjointness of the Schrodinger operator is now ensured 
by imposing the unitarity of the scattering matrix : Q^Q = 1. The relation between the two 
formulations yields : 

Q = (^D-C)- 1 (^D + C). (54) 

For the discussion below it is convenient to introduce the parameter 7 related to the energy by 
E = —7 = k 2 + i0 + (the matrix Q is unitary for E > only) . An example of a matrix Q for one 
vertex without continuity of the wave function was given in refs. |1921 HUtfl I52j . This particular 
choice has become popular in mesoscopic physics (see also ref. [209! where a more convenient 
parametrization was provided). 

Since we are considering here compact graphs, we are dealing only with discrete spectrum. The 
study of noncompact graphs (with some wires of infinite length) , which have continuous spectra, 
requires a scattering theory approach, initiated by the work of Shapiro (1911 1192j and discussed 
in refs. fTTTRl 1^771 01 rTTTT?! fTl^l I^Kl fTim I^THTl I^TF^l l^nni 077^] . 

Magnetic fluxes.— The Laplace operator arises in the context of diffusion equation but also in 
quantum mechanics. In this last case a natural generalization is to introduce a magnetic field. 
This is achieved by introducing a 1-form A(x)dx along the wires (the derivative must then be 
replaced by a covariant derivative : d x — > D x = d x — iA(x)). We denote by 9 a p = f^dxA(x) 
the corresponding line integral along the arc af3. In the context of classical diffusion, magnetic 
fluxes and winding numbers are conjugated variables. 
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5.3 Trace formulae and zeta functions 



Trace formulae play an important role in spectral theory. A famous example is the Selberg trace 
formula which may be viewed as an extension of the Poisson summation formula to non commu- 
tative groups |190| . An analogous formula in physics is the Gutzwiller trace formula 123^ that 
has been extensively used in the context of quantum chaos and mesoscopic physics. Although 
one is exact and the other only a semiclassical approximation, both of them express the partition 
function (or the density of states, whose Laplace transform gives the partition function) as a 
sum over closed geodesies. They provide a connection between quantum properties (spectrum) 
and classical properties (classical trajectories). Below we discuss two examples of exact trace 
formulae that have been derived for graphs and their relation to spectral determinants. 



Roth's trace formula 

It expresses the trace of the heat kernel (partition function Z(t) = Tr{e* A }) as an infinite 
series of contributions of periodic orbits on the graph. This remarkable formula, due to Roth 
|186l 1185] . applies to graphs with continuous boundary conditions with X a = 0. It it easy to 
include magnetic fluxes additionally : 

Z(t) = ^r= + ^ + ^r=X KC)a(C)e-^^ . (55) 

ZyTTt ^ zV Tit 

C is the "volume" of the graph, i.e. the total length C = Ylfa^^ap- The sum runs over an 
orbits C = (h,t2, • • • ,in) constructed in the graph. 1(C) = + • ■ ■ + k n is the total length of 
the orbit, and 6(C) the magnetic flux enclosed by it. C designates the primitive orbit associated 
with a given orbit C. The weight a(C) depends on the connectivity of the vertices visited by the 
orbit : a(C) = e^e^ • • • The matrix e couples the arcs of the graph 30 : 

• if i ends at vertex a and j starts from it, we have = 2/m a , where ra a is the connectivities 
of the vertex. 

• If moreover i and j = i are the reversed arcs eg = 2/m a — 1. 

• Otherwise = 0. 

It is worth mentioning that the Roth trace formula has found recently some pratical appli- 
cations to analyze magnetoconductance measurements on large square networks 31 94 . 



Ihara-Bass trace formula 

Instead of considering metric graphs we now turn to graphs viewed as purely combinatorial 
structures consisting of vertices connected by bonds of equal lengths (l a p = 1). All the infor- 
mation is therefore encoded in the adjacency matrix A. In this setting, the Ihara ^-function is 
defined as 

cw- i =n( i -^ (c ~ fl) ) > (se) 

C B 

where the infinite product extends over all primitive backtrackless orbits Cb- The Ihara-Bass 
trace formula relates this infinite product to the determinant of a finite size matrix [281 1197j : 

((u)" 1 = (l - u 2 ) B ~ V det((l - u 2 )l -uA + u 2 Y) (57) 

30 Note that the parameter has a simple interpretation in the scattering formulation |2U9| : it is the probability 
amplitude to be transmitted from arc i to arc j. The arc matrices e and Q are related by = where j 
denotes the reversed arc. 

31 More precisely, the expansion given below by eq. (1621 for the Laplace transform of the partition function is 
relevant in this case. 
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where 1 is the identity matrix, A is the adjacency matrix (A a p = a a p) and Y is the diagonal 
matrix encoding all connectivities : Y a R = 5 a pm a . This relation was derived by Ihara |129j for 
regular graphs (all vertices with same connectivity) and later generalized to arbitrary graphs by 
Bass |28| . The formalism that we have developed for metric graphs is flexible enough to describe 
these combinatorial structures in the same setting. We will see below that the ^-function is in 
fact directly related to the spectral determinant. Moreover the Ihara-Bass trace formula can be 
further generalized to include backtrackings. 

5.4 Spectral determinant 

In the physics literature spectral determinants arise in evaluating path integrals which are 
quadratic in the fluctuation around a given background field. A well known technique for 
regularizing such quadratic path integrals is the ^-function regularization. Given a certain oper- 
ator O whose eigenvalues E n are known, one defines the following ^-function : £(s) = ^2 n E~ s . 
This expression, which converges for s sufficiently large, can be analytically extended to a mero- 
morphic function regular at the origin. The corresponding regularized determinant is then 
det O = exp(— C'(0)) (= Fln^n formally). One can find a general discussion on functional 
determinants in ref. |96| . 

In the context of graphs another regularization of the determinant of the Laplace operator 
has been used and has proved to be directly related to several physical quantities. If we introduce 
the trace of the resolvent 5(7) = ^ n (7 + En) -1 , the spectral determinant is defined as ^(7) = 
exp(J 7 d7 ; 5(7')) (it can be formally written as 32 5(7) = Iln(7 + E n ))- Moreover, the spectral 
parameter 7 has in some cases a physical meaning (see section l5~6|) . 

Laplace operator A with continuous boundary conditions 

Pascaud & Montambaux have shown in refs. |177| 1178] that S(-y) = det (7 — A) can be related 
to the determinant of a V x ^/-matrix : 



where the product runs over all bonds of the network. The V x ^-matrix M is defined as : 



where the adjacency matrix constrains the sum to run over all vertices /i connected to a. The 
matrix M encodes all information about the network : topology (matrix a a p), lengths of the 
wires (l a p), magnetic fluxes (6 a p), boundary conditions (X a ). Below, we give several examples 
which show how S(j) is related to the characteristic function of various interesting functionals 
of Brownian curves on a graph. 

The expression (|58[1 was originally derived in ref. |178j by constructing the Green's function in 
the graph and eventually integrating it : f dx ( x | \x) = In S(j) . A more direct derivation 
using path integral was later obtained in ref. jHJ. 

The relation between the Roth trace formula and the result (|58j) was addressed in ref. [Hj. The 
main difficulty to establish this connection is to go from vertex variables to the arc language 
of eq. (|55|). A first step is to relate the determinant of the vertex-matrix, eq. (|58|) . to the 
determinant of an arc-matrix : 



S(i) = 7^ 11 sinh(^/ a/3 ) det M , 

(a/3) 



(58) 




(59) 



5(7) = 7— e^ c det(l - QR) . 



(60) 
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The limit 7 — > of S(y) has been discussed in ref. 0- 
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R is the 2B x 2B matrix Rij = 5^ e~^ li+l9 ' , where j denotes the reversed arc. The matrix Q 
was introduced above (eq. (|54|l) and is related to e by Qij = e^. For the continuous boundary 
conditions with \ a = it is given by Qa = 2/m a — 1, Qij = 2/m a if i and j both issue from 
the vertex a. Qij = in other cases. Eq. H6U|) holds for the most simple boundary conditions : 
continuous with X a = 0. The general case is discussed below. Expanding the determinant by 
using lndet(l — QR) = — YH?=i h Tr {((J-R)" - }, we eventually express the spectral determinant 
as an infinite product over the primitive orbits : 

5( 7 ) = ^e^ c H (l - a(C)e-^ +w 0) . (61) 
c 

This shows that the spectral determinant is a zeta function (references on zeta functions on 
graphs are refs. |197| lB7|). The last step to connect this formula to Roth's trace formula is to 
notice that 

- ir, + ^EKcnne-^-m , (62) 

where the sum now runs over all orbits (if C is not primitive, C designates the related primitive 
orbit). Finally we perform an inverse Laplace transform of this expression, dt Z(t)e~~ /t = 
-Tp-lnS (7) and eventually recover eq. ((531) . Examples of applications are studied in ref . j3] . 

Schrodinger operator —A + V(x) with general boundary conditions 

The result (|58|) of Pascaud & Montambaux has been generalized by one of us. In refs. [721 173j a 
similar formula was obtained for the spectral determinant of the Schrodinger operator — A+V(x) 
(the Hill operator), where V(x) is a scalar potential defined on the graph. In ref. [71] the 
formula was further extended to describe general boundary conditions as well. As an illustration 
we construct the generating function of the number of closed orbits with a given number of 
backtrackings (1231 ISSH] . 

The starting point is to introduce two linearly independent solutions of the differential equa- 
tion (— d^ + V a p(x) + j)f(x) = on [0, l a p\. We associate each solution with an arc. Let us 
denote f a p{x a p) the function satisfying 

f af3 (0) = l and f a p(l a p) = 0. (63) 

Therefore a second solution of the differential equation is naturally denoted f/3 a (%f3a) = f/3a(l>a/3 — 
x a! 3). The Wronskian of these two solutions, defined as W a p = / a ^(i a g) d fp^ x P^ _ A fap{x afi ) ^ a ^ Xj3c 
is constant along the bond : W a p = Wp a = —f' a g(lap) = — /g a (W)- ^ we consider the case 
V(x) = 0, the solution is simply f a ? M = - S^g- 

All the required information about the potential is contained in the 2B x 2B arc-matrix N, 
defined as 

NaPM = S afJ ,S/3 v f' aj3 (0) - 8 arj 5p^f al3 {l a i3) . (64) 

This matrix couples a given arc to itself and to its reversed arc, only. If we assume that the 
matrices C and D are independent on the spectral parameter 7, it was shown in ref. |74| that : 

5(7) = det( 7 - A + V{x)) = TT det(C + DN) (65) 

(a/3) WaP 

where the product runs over all bonds. Functional determinants on a segment of R with general 
boundary conditions at the boundaries have been studied by McKane & Tarlie j!64j using the 
formalism developed by For man |96j . 
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as 33 : 



It is also interesting to encode the information on the potential V{x) in the matrix R defined 

R={^l + N){^jl-N)- 1 . (66) 



Then 

s « = n 4^ ^drw) det(c - det(1 - QR) (67) 

where the matrix Q = (y/^fD — C)~ 1 (t [ /jD + C) was defined above by eq. (|54|). Let us explain 
the structure of eq. ((57)1 : the term det(l + it!) II(a/3) contains informations on the potential 
only 34 . The factor det(C — 1/7D) contains only information on the topology of the graph. The 
most interesting part is the last term det(l — QR) combining both informations. In particular 
this last part generates the infinite contributions of primitive orbits in 1)6 If) • 

Permutation-invariant boundary conditions. 

In the previous section the matrices describing the graph are 2B x 2B arc matrices. A simplifi- 
action can be brought by passing to vertex variables. In the arc formulation, C and D define the 
topology of the graph : two arcs coupled by C and /or D issue from the same vertex. Therefore 
it is possible to organize the basis of arcs in such a way that the matrices C and D have similar 
block diagonal structures. The matrices are made of V square blocks, each corresponding to a 
vertex. A given block, of dimension m a x m a and denoted C a (and D a ), corresponds to the m a 
arcs issuing from the vertex a. If we assume that the boundary conditions are invariant under 
any permutation of the nearest neighbours of a, then it is possible to introduce vertex variables. 
In this case we can write: 

C a = c a l + t a F a (68) 
D a = d a l + w a F a (69) 

where F a is a matrix with all its elements equal to 1. The boundary conditions at the vertex 
a are characterized by the four parameters c a , d a , t a and w a (note however that this choice is 
not unique). 

Now, let us show that, for boundary conditions given by eqs. I)68I69J) and V(x) / 0, the 
spectral determinant can be expressed in terms of the vertex V x ^-matrix. 

We proceed as before but, this time, we consider, for each bond, two other independent 
solutions, Xa/3(x a /3) and X/3a(l a f3 - x a p) = X/3a(xp a ), of the equation (-d| a(3 + V a p{x a p) + 
7) x{ x aB) = that satisfy the following conditions: 

c a x*p(0) + d a x' af i(0) = 1 (70) 
CpXapilap)- dpXapttap) = (71) 

We denote by H a/ g the Wronskian of Xa/3 and Xfia- Following the same steps as before, we get 
the spectral determinant (up to a multiplicative constant) : 

5(7) = Yl J— det M (72) 



33 Let us remark that, for the free case (V(x) = 0), we recover the expression of the matrix R given above, 
coupling the arc af3 to the reversed arc /3a only : Rap,^ = SarjSp^e - ^ 10113 . 

34 The matrix iV encodes the information about the potential on the bonds through /4a(0) and f' a p(lap)- This 
information can also be introduced through transmission t a @ and reflection r a p amplitudes by the potential 
V a p(x). This relation is developed in ref. |209| . It is interesting to point out that det(l + R)Yi( a f3) W a p = 

n B B/2o 4. oS S /2 tt D 

2 7 li( Q/ 3)W>=2 7 ll (Q/3) -R Q /3,/3 a . 
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where M is the V x V-matrix : 

M a/3 = 5 aj3 ^1 + ^ a atx [t a Xan(0) + w a x' afl (0)} J + a al3 [c a w a - t a d a ] E a/3 . (73) 

The eq. 1)65(1 . expressing the determinant in terms of arc matrices, and the eq. ()72j). expressing 
it in terms of vertex matrix, have been derived up to multiplicative constants independent on 7. 
We can establish a precise relation by comparing their behaviour for 7 — > 00. We end up with : 

[] det(C + ZW) = [J J- det M . (74) 

Free case (V(x) = 0). Applications : counting backtrackings 

In this subsection we show an application of generalized boundary conditions to count back- 
trackings (the figure shows an orbit with one backtracking). 

We now consider the case V{x) = still with permutation-invariant boundary conditions. 
With the notations 

eqs. (HI ESI El lead to : 

det(l - QR) = 2- y H (p aVa ) H(l- ^e-VrW) d et M (75) 

a (a/3) 

and the V x ^-matrix M takes the form : 

/ 2 m„ 1 ^ 1 + 2e-^ 

Note that the expression (|51?|) is recovered for r\ a = \ and = 2/(m a + X a /yfy). 

For permutation-invariant boundary conditions, the matrices C, D and Q (eq. 1)54(1) are 
block-diagonal. The block Q a takes the simple form: 

Qa =Va(-l + PaF a ) (77) 

The only non- vanishing elements of the matrix QR are : 

(QR)a(3, m = {PaVa ~ Va Sfr) e~^ la » (78) 



We call paJ]a the transmission factor at vertex a. and p a i]a ~ Va the reflection factor. By using 
the same expansion as above we write 

det(l -QR)=Y[(l- p(C)e~^) , (79) 
c 

where the product is taken over all primitive orbits C whose lengths are denoted by /(C). An 
orbit being a succession of arcs . . . , ra, aft, . . . with, in a, a reflection (if r = (3) or a transmission 
(if t ^ (3), the weight //(C), in eq. (|79j). will be the product of all the reflection - or transmission 
- factors along C. 
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Graphs with wires of equal lengths. - We consider the case of equal lengths l a p = I ; then we can 
choose 7 = 1 without loss of generality and introduce the notation u = e~ l . It is clear from 
eq. (|77|) that a backtracking at vertex a brings a factor r/ a in the weight /u(C). In order to count 
backtrackings one has to choose the boundary conditions p a r) a = 1 and rj a = rj. Eq. ((751) takes 
the simple form : 

Yl (l - (1 - r]) nR ^u m ^ = (1 - ifu 2 ) B - y det ((1 - f?V)l + r]u 2 Y — u A) = Z' 1 (80) 

Cm 

m is the number of arcs of the primitive orbit C m and nn(C m ) is the number of reflections 
(backtrackings) occuring along C m . Y is the V x F-matrix Y a p = 5 a p m a and A is the adjacency 
matrix. 

Setting r] = 1 implies riR(C m ) = in the left-hand side of eq. (|80j) : we recover Ihara-Bass 
formula 129, 1281 1197j where only primitive orbits without tails and backtrackings are kept. 
(Ihara |129j established this formula for a regular graph ; the proof for a general graph is given 
in refs. |28l I197j using a direct counting technique). 

Now, let us consider closed random walks with a given number of backtrackings. 

Eq. (|8Ujl provides a non-trivial generalization of the Ihara-Bass formula (an independent 
derivation is also given in ref . 26 ) . As an application let us consider the problem of enumerating 
m-steps random walks with p-backtracking steps |226j . Taking Z in (|80|). we get: 

j , y oo m V 

u ^Jr = E E E Ar ™(«) c 1 - v) p u m (si) 

m=2 p=0 a=l 

where Nm(a) is the number of m-steps closed random walks on the graph starting at a, with p 
backtrackings. 

Example : The complete graph. - For the complete graph 35 Ky, we get the results: 

l)(V-2) 
l)(V-2)(V-3) 
1)(V - 2)(V - 3)(V - 4) 

l)(y-2)(T/ 3 -9F 2 + 29V-32) (82) 



]\£(a) = Nl(a) = 
5(V-l)(V-2)(V-3) 

6(y-i)(y-2)(y-3) 2 (83) 

Note that these expressions have been obtained when vertices are all characterized by the same 
parameter rj. Introducing different parameters r\ a allows counting the backtrackings at a given 
vertex. 

35 The complete graph Kv with V vertices is the V — 1-simplex : each vertex is connected to all other vertices 
(see figure 



iv 2 V) = o 
iv 3 V) = (v- 

K(a) = (V — 
N$(a) = (V- 
N°(a) = (V- 

and also: 

Nfra) = 

N£(a) = 

N*(a) = 
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Figure 6: The complete graph K§ (the 4- simplex) : all vertices are connected to each other by 
bonds of equal lengths. 



5.5 Quantum chaos on metric graphs 

It has been recently realized that metric graphs are interesting models for quantum chaos. The 
paper of Kottos & Smilansky (140] has stimulated several works on spectral statistics and level 
correlations [1411 1241 12L)1| , and renewed the interest in the Roth trace formula |186j . Progress 
in the understanding of universality of spectral statistics for generic quantum graphs has been 
achieved by Gnutzmann & Altland |11U| . One of the input of this work is the observation that 
the spectral average for a given graph with incommensurate bond lengths is equivalent to an 
average over a certain ensemble of unitarity matrices. Star graphs |36U136j have provided simple 
examples of systems with intermediate level statistics similar to the one observed in Seba billiards 
and a precise connection has been established in ref . . There is however an interesting class of 
graphs which do not enter in this category |2U2| . Chaotic scattering and transport properties in 
open graphs (with some infinitely long wires) have also been studied in refs. |142[ I25 |I143| . Since 
quantum chaos is not the central subject of our review, this list of papers is not exhaustive : 
we refer the reader interested in this topic to the review papers |141( I143| I202j and the PhD 
thesis [22S] , 



5.6 Weak localization and Brownian motion on graphs 

Our initial physical motivation for the study of the spectral determinant on graphs was based 
on the observation that the weak localization correction to the conductivity is directly expressed 
in terms of the spectral determinant. We see from eq. (j2J that 

(Aa} = - 2 f ^ In S( 7 ), (84) 
7T Volume c*7 

where the spectral parameter is related to the phase coherence length 7 = 1/ 1? . The expressions 
(|84I58I59|) , due to Pascaud & Montambaux [1771 11781 |3] , improve the approach initiated by 
Dougot & Rammal |HT1IH2] . 



Nonlocality of the quantum transport in arbitrary networks 

We must stress that the formula (|84jl corresponds to a uniform integration of the cooperon, 
defined as P c (x,x) = (x\-^\x), on the graph (Act) oc — f dx P c (x,x). This approach is 
limited to the case of regular graphs, where all wires play the same role. In other terms, 
for a nonregular network, the quantity (|84[) does not correspond to a quantity measured in a 
transport experiment. For arbitrary networks, the cooperon must be integrated over the network 
with appropriate nontrivial weights that depend on the topology of the whole network and the 
way it is connected to external contacts. Such a generalization was provided in ref. |210| . 
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Windings in a loop connected to a network 

As we have mentioned in the section 12.21 magnetoconductance oscillations due to a magnetic 
flux and winding properties are closely related. For example, if we consider an isolated ring of 
perimeter L pierced by a magnetic flux </>, the well-known behaviour of the harmonics 

(Aa n )= / -{Aa(e)}e- rafl oce- |n|L ^ (85) 
Jo 27r 

is a direct consequence of the fact that the winding around the ring scales with time as nt oc 

t l/2 

(normal diffusion). 9 = 4ir^>/<po is the reduced flux and 7 = l/L^. This effect was predicted 
by Al'tshuler, Aronov & Spivak (AAS) in ref. ^U] and observed in experiments on cylinder 
films [1331 ITT]. 

Recently it has been noticed that the fact that the ring is connected to arms, which is 
necessary to perform a transport experiment, can strongly affect the harmonics. If we consider 
a ring connected to N a long arms, in the limit L v <C L, the harmonics are still given by the 
AAS behaviour (|85[). however when the perimeter is smaller than L„, the harmonics behave 

as (Aa n ) oc e~\ n \V NaL / L v . To clarify the origin of this behaviour, the winding of Brownian 
trajectories around a ring connected to another network has been recently examined in ref. [212] . 
This analysis is based on the fact that the winding number distribution can be expressed in terms 
of the spectral determinant. Let us briefly describe this approach. 




Figure 7: A diffusive ring attached to a long arm. The parameter Ai describes the boundary 
condition at vertex 1 : X\ = 00 for the Dirichlet boundary and X\ = for the Neumann boundary 
condition. 




Figure 8: A diffusive ring attached to (a) : an infinite square lattice, (b) : an infinite Bethe 
lattice of connectivity z = 3. To ensure that the Bethe lattice is regular, it must be embedded in 
a constant negative curvature surface. 

We consider a ring to which is attached an arbitrary network at vertex (the figures and 
[H]give examples of such a situation). Our aim is to understand how the winding around the ring 
is affected by the presence of the network. To answer this question we introduce the probability 
to start from a point x on the graph and come back to it after a time t, conditioned to wind n 
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times around the loop : 

rx(t)=x 

V n (x,t\x,0)= Vx(T)e-y° dT± ^ 2 5 nM[x{T)] (86) 

Jx(0)=x 

where Af[x(r)] is the winding number of the path x(t) around the loop. For simplicity, we 
consider the case where the initial point is the vertex 0. The computation of "P n (0, t\0, 0) requires 
some local information (the eigenfunctions of the Laplace operator at point x = 0). On the 
other hand the spectral determinant encodes a global information since it results from a spatial 
integration of the Green's function of the Laplace operator over the network. However we have 
shown in ref. |212j that if we consider mixed boundary conditions at vertex described by a 
parameter Ao (see the definition of the boundary conditions given above), we can extract the 
probability of interest 36 as follows : 

POO p2ix J 

/ dtV n (0,t\0,0)e-^= / e~ ine —InS^ij) , (87) 
Jo Jo 27r dA a =o 

where S^°\j) is the spectral determinant for mixed boundary condition at vertex and with 
a magnetic flux 6 piercing the ring. This formula allows us to express the probability as 

I 00 df P B (0,t|0,0) e"* = J- . f h pf L , , e-K^-M , (88) 
Jo 2 v Aysmh v AyL eff (7) 

where all the information about the nature of the network attached to the ring is contained in 
the effective perimeter L e g (7), defined as : 

Slim a/7-L 

cosh VT^eff (7) = cosh yfyL + — — T . (89) 

1 l 7W netJ00 

The matrix M ne t describes the network in the absence of the loop (it is given by eq. (|59|) ; for the 
network of figure |EH> A^net describes the infinite square network without the ring). The matrix 
element (M~J^)qq has a clear meaning : it is the Green function of the Laplace operator in the 
network (without the loop) computed at the position where the ring is attached : v / 7(Af r ^,J;)oo = 

(QI^aIo). 

The effective perimeter L e g{~f) probes the winding at time scale t ~ I/7 : precisely, the 
winding number scales with time as nt ~ Vi/L e g(l/t). 

Two interesting examples are 

• Ring attached to an infinite wire (figure^.- When a wire of length b with the Dirichlet 
boundary at one end is attached to the ring, it is easy to see that eq. ((S9*|) gives (M^)qo = 
tanh^/7?) — ► 1. At large time t 3> L 2 , the effective perimeter behaves as L e ff ~ \/^7 _1//4 . This 
behaviour is related to a scaling of the winding number with time 

n t oc i 1/4 . (90) 

The full distribution for the winding number n is given in ref. |212| . The exponent 1/4 that 
characterizes anomalously slow winding around the loop originates from the fact that the dif- 
fusive trajectory spends a long time in the infinite wire, which increases the effective perimeter 

,j6 The fact that some "local" information, such as "Pn(0, t|0, 0), can be extracted from a more "global" object, 
like the spectral determinant, by introducing the mixed boundary conditions with parameter Ao has been also 
used in the context of scattering theory in graphs in refs. |2()rill2()^] . The derivative in eq. I|87|l can be understood 
as a functional derivative since Ao plays a role similar to the weight of a 5-potential at xq. The use of functional 
derivatives in scattering theory for mesoscopic systems has been fruitfully used by Biittiker and co-workers (see 
ref. |51) and references therein). 
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at such time scales. This problem is also related to the anomalous diffusion along the skeleton 
of a comb, studied in refs. [2221 122j by different methods. It is interesting to use this picture. 
Let us consider a random walk along the sites of a line where the diffusive particle is trapped 
during a time r on each site. The trapping time is distributed according to a broad distribution 
Pi(t) oc r ^ with < /J, < 1. It follows that the distance scales with time as n t ~ W 2 02]. If 
we go back to the problem of diffusion along the skeleton of a comb (or the winding in the ring 
connected to the long arm), the arm plays the role of the trap. The distribution of the trapping 
time is given by the first return probability of the one-dimensional diffusion : P\(r) oc t~ 3 / 2 
and we recover eq. (|90|) . 

• Ring attached to a square network (figure\$p,).- When studying the winding around the loop, 
it is important to know whether the Brownian motion inside the network attached to the ring 
is recurrent or not. Let us consider the case where the network attached is a d-dimensional 
hypercubic network. For d > 2 the Brownian motion is known to be transient where as, for 
d = 2 it is neighbourhood recurrent (in d = 1 the Brownian motion is pointwise recurrent). 
Therefore we expect the dimension 2 to play a special role. In the large time limit, when t S> L 2 
and t S> a 2 (a is the lattice spacing), we find an effective length ^/j L c g ~ \J ain(4/^/ ya)> *hat 
corresponds to a scaling of the winding around the loop 

n t oc (lnt) 1 / 2 . (91) 

This result can be obtained in the same way as for the ring connected to the arm. This time the 
plane acts as a trap. The distribution of the trapping time is given by the first return probability 
on a square lattice, which is known to behave at large times like P\(r) oc l/(rln 2 r) [211, from 
which we can recover eq. (|91j). 

Occupation time and local time distribution on a graph 

Another set of problems concern occupation times : i.e. time spent by a Brownian particle 
in a given region. An example of such a problem is provided by the famous arc-sine law for 
the Id Brownian motion that gives the distribution of the time spent by a Brownian motion 
(x(t),0 ^ r ^ 1 1 x(0) = 0) on the half line M + . This result was derived long ago by P. Levy 
|155j . It has been extended by Barlow, Pitman & Yor [23] for a particular graph (star graph 
with arms of infinite lengths) : instead of an infinite line, these authors consider n semi-infinite 
lines originating from the same point and study the joint distribution of the times spent on each 
branch. More recently, this problem has been reconsidered in the case of arbitrary graphs 75 . 
It may be stated as follows : consider a Brownian motion x(r) on a graph, starting from a point 
xo at time and arriving at a point x\ at time t. Let T a n denotes the time spent on the wire 
(a(3). This functional is defined as T a p[x(r)] = L dr Y a p(x(T)) where the function Y a p(x) is 
1 for x S (cxf3) and otherwise. Our aim is to compute the Laplace transforms of the joint 
distribution 

(e-EwtapT^ = J tx^ixuxo-MZap}) (92) 

with 

rx(t)=Xi t 

F(x 1 ,x ;t;{£ a p})= Vx(t) e"^ dr± . (93) 

J x(0)=xo 

The definition (|92jl takes into account averaging over the final point x\. A conjugate parameter 
is introduced for each variable T a p, i.e. each bond. A closed expression of the Laplace 
transform of the joint distribution (|92|) has been derived in ref. [73]. The result is given as a 
ratio of two determinants, an expression reminiscent of the one of Leuridan |154j . although the 
connection is not completely clear. 
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Similar methods have been also applied in ref. jSH] to study the distribution of the local time 
T x [x(t)] = J *dr5(x(r) - x ). 

A simplification occurs when the final point coincides with the initial point. Then it is 
possible to relate the characteristic function to a single spectral determinant. We do not develop 
the general theory here but instead consider an example close to the one studied in the previous 
subsection. Let us consider the graph in figureEland ask the following question : for a Brownian 
motion starting from xq at time and coming back to it at time t, what is the distribution of 
the time T' arm [a;(r)] spent in the arm if in addition the Brownian motion is constrained to turn 
n times around the ring ? It is natural to introduce the following function : 

px(t)=xo t 

F n (x ,x ;t,0= / Vx(T)e-fo*TU***™(*)) 6 nMx{T)] (94) 

J x(0)=xo 

where JV[:e(t)] is the winding number around the ring. This function is related to the Laplace 
transform of the distribution of the functional T arm [x(r)] 

-fT«m[x]\ _ Fn(xQ,X ;t,£) ^ 

I c n F n (xo,x Q ;t,Q>) ' 

where (• • -) c denotes averaging over curves of winding n. The denominator ensures normaliza- 
tion. The Laplace transform of (|94|) is given by a relation similar to eq. (|87|) 

f°° C 2lr r\f) H 

/ dtF n (x ,x ;t,Oe-^= e-^-^ln^')( 7 ) (96) 

Jo J 2vr dA A =o 

where the appropriate spectral determinant is built as follows : (i) since the starting point is 
fixed at 0, we introduce mixed boundary conditions at this point, with a parameter Ao that will 
be used to extract the "local information" . (ii) A magnetic flux 9 is introduced (conjugate to the 
winding number). (Hi) The spectral parameter is shifted in the arm as 7 — > 7 + £ to introduce 
the variable £ conjugate to the time T arm . 

We choose the vertex as initial condition and impose the Dirichlet boundary condition 
(which is achieved by setting Ai = 00) at the end of the arm of length b. The spectral determinant 
is found straightforwardly (for an efficient calculation of S(j) for a graph with loops, see ref. |177j 
or appendix C of ref. 0) : 



S(*>)( 7 ) = fi nh ^ L sinhV7 + f & 



cosh \/7-L — cos 8 , , 

2 ^ sinh^yL + A + \/7 + £coth y^ + f & 



(97) 



The terms in the brackets correspond to the matrix 3 M. (since Ai = 00 we can consider only 
the element .Moo)- The first term is the contribution of the loop, the second comes from the 
boundary condition and the last one comes from the arm. It immediately follows that 

^ dt F n (0, 0; t, e-* = ^- Si ff T L e-W^ (98) 
Jo 2 A /7Sinh A /7L eff 

with ^ 

cosh ^L cS = cosh ^jL + - yT+ £/ 7 sinh coth a/7+~£ b . (99) 

Let us consider an infinitely long arm b — > 00. If we are interested on time scales t 3> tl, 
where tl = L 2 is the Thouless time over which the ring is explored, eq. (|99|) gives ^/7-L e fj ~ 
(7 + £) 1/4 \/Z, therefore from eq. © we see that ^"„(0 5 0;t, - ^ ^n(0, 0; t, 0). The inverse 
Laplace transform of eq. (|95|) leads to (S(T — r a rm[^]))c n — S(T — t), which means that the 
Brownian motion spends almost all the time in the arm (this simple result confirms the picture 
presented to explain the scaling of the winding around the ring of the form nt oc i 1//4 ). 

37 The introduction of the conjugate parameters {£0/3} corresponds to shift the spectral parameter 7 on each 
wire as 7 — > 7 + £ a p. The matrix to be generalized is not M, given by eq. JHHJ, but M = y/yM. 
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6 Planar Brownian motion and charged particle in random mag- 
netic field 



A model of random magnetic field describing a charged particle moving in a plane and subjected 
to the random magnetic field B{r) = V x A{r) due to an ensemble of magnetic Aharonov-Bohm 
vortices is described by the following Hamiltonian |77| 1100] : 

H = \ (jf- Mrl) 2 + \ B{f) = \ [v- «E } ) +^J25(r-r i ) (100) 

where ^ a < 1. We have set H = e = m = 1. a is the magnetic flux per vortex, in unit of the 
quantum flux 0o = h/e : a = (f>/(f>o = 4>/(2ir). The model is periodic in a with a period 1. u z 
is the unit vector perpendicular to the plane. r*j are the positions of the vortices in the plane : 
they are uncorrelated random variables (Poisson distribution). The 5 interactions (coupling to 
the magnetic field) are necessary in order to define properly the model |131[ I163[ I175| 133"! 163] . 

By comparison with scalar impurities discussed in the introduction, the magnetic nature 
of the scatterers gives rise to rather different properties 38 : it was shown in ref. (23 that the 
spectrum is reminiscent of a Landau spectrum with Landau levels broadened by disorder, in 
the limit of vanishing magnetic flux a — ► (see below). This analysis was performed using 
perturbative arguments supported by numerical simulations. The latter use a relation between 
the average DoS of this particular model and some winding properties of the Brownian motion. 

Let Z(t) denotes the partition function for a given distribution of vortices and Zq{£) the 
partition function without fluxes. The ratio of partition functions can be written as a ratio of 
two path integrals 



Q l fc A ' ar \ (101) 

where Zq{€) = V/(2nt) is the free partition function (V is the (infinite) area of the plane) and 
(■ ■ - )c stands for averaging over all closed Brownian curves of the plane. 

In order to average Z(t) over the Poissonian distribution of vortices, let us consider, for a 
while, our problem on a square lattice with lattice spacing a. Let Ni be the number of vortices 
in square i. the magnetic flux through any closed random walk C on this lattice can be written 
as : 

A - &r = ^ j 2iraN i n i (102) 



C 



where rii is the number of times the square i has been wound around by C. 

Averaging Z(t) with the Poisson distribution (// being the mean density of vortices) 



(/ia 



2\Ni 



gives 



PW) = ^-e-" a (103) 



(Z(jt)) w = Z (t) (e^ Xdf? ) c{ ., } = Z Q (t) / exp L^S n (e 2i ^ - 1) ) ) (104) 




38 For scalar impurities, the presence of a magnetic field also strongly affects the spectral properties. The case 
of a Gaussian disorder projected in the lowest Landau level (LLL) of a strong magnetic field has been studied in 
ref. |220| . Other disordered potentials have been considered later. In particular it was shown in ref. 46i that for 
a weak density of impurities, the spectrum can display power law singularities. A physical interpretation of such 
power law singularities has been discussed by Furtlehner in ref. |100l HUl) . The relation between the model of 
scalar impurities projected in the LLL of a strong magnetic field and the model of magnetic vortices was discussed 
in refs. |78l I1UU| . Finally we mention a recent work on Lifshitz tails in the presence of magnetic field |153| . 
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where (• • -){rj denotes averaging over the positions {r*j} of the vortices. 

The quantity S n denotes the area of the locus of points around which the curve C has wound 
n times. This result was derived for a random walk on a lattice but is also obviously valid 
off lattice as well, for Brownian curves. The winding sectors of a closed curve are displayed in 
figure El 








Figure 9: A closed curve with its n-winding sectors. The label n of a sector is the winding 
number of any point inside this sector. 



The random variable S n scales like t. In particular, we know |62U118| its expectation (S n ) = 2 ^ . 
Moreover, it has been shown in 223 that the variable n 2 S n becomes more and more peaked 
when n grows: P(n 2 S n = X) — > 5(X — (the average area of the n = 0-sector has been 

n— >oo 

recently studied in ref. jl(J4j ). 

Thus, extracting t, eq. (|1U4|) rewrites as 

(Z(t)) m = Z (t) J dSdAP(S,A) e-^ 5+iA ) = Z (t) (e~^ s+iA ^ c , (105) 

where P(S, A) is the joint distribution of the rescaled (t independent) variables S and A defined 



as 



2 

S = - S n sin 2 (-7ran) (106) 

n 

A = - S n sin(27ran) . (107) 



t 

n 

The averages of these two variables are given by : 

(S) = vra(l-a) (108) 
(A) = 0. (109) 

With eq. ()105|) . we observe that y (Z(t))| r -?j has the scaling form F([it)/t. Thus, its inverse 
Laplace transform, the average density of states per unit area, is a function of only E//J, and a. 
Moreover, it is easy to realize that (#(i))r-i is even in a (each Brownian curve in {C} comes 
with its time reversed) and periodic in a with period 1. Thus, one can restrict to ^ a ^ I. 
Let us focus on the two limiting cases a — > and a = 1/2. 

• Limit a — ► : Landau spectrum. 

When a — > 0, a careful analysis shows that 

(Z(t)) m * Z {B} e-iW* (110) 
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where = ^o(0 S i n h((g)t/2) * s ^ e Petition function of a charged particle in a uniform magnetic 
field (B) = 2npa (we recall that the impurity i carries a magnetic field 2ira5{r — r*j)). The system 
of random vortices is, thus, equivalent to the uniform average magnetic field, albeit with an 
additional positive shift in the Landau spectrum : the inverse Laplace transform of the partition 
function <fTT7^ gives the Landau spectrum made of equally spaced infinitely degenerated levels. 
This corresponds to the oscillating behaviour shown on figure 1111 The origin of the shift can 
be traced back to the presence of the repulsive 5 interactions that have been added to the 
Hamiltonian to define properly the model. 

• Half quantum flux vortices (a = 1/2) — The spectral singularity at E = 0. 

In this case the variable (|1U7|) vanishes, ^4 = 0, implying that {Z(t))r^,x = Zo(t){exp (—p,tS))c 
where now 

S=*^S n . (Ill) 

n odd 

Performing the inverse Laplace transform, we get the average density of states per unit area: 

(p(E))=p (E) / dSP(S) (112) 
J o 

where po(E) = ^- is the free density of states per unit area and P(S) is the probability distri- 
bution of S. 




Figure 10: The average density of states at a = 0.5 (the free density of states is constant 
Po = 1/(2^)^). Circles are simulation results. The full line is the fit discussed in the text. 

We may use this result to determine the nature of the singularity in the average DoS. As 
displayed in figure ITU1 (where we have taken p = 1), {p(E)) increases monotonically from 
to po(E) with a depletion of states at the bottom of the spectrum. Circles are the result of 
numerical simulations (on a 2D square lattice, we have generated 10000 closed random walks of 
100000 steps each). The full line is a low-energy fit of the quantity (p(E))/po by the function 

P^l= ae -K,/Ef (113) 

PO 

with a = 2.8 and b = 1.4 (see figure HU|) . It is worth stressing that this behaviour is quite 
different from the one expected for a disorder due to a scalar potential. The famous Lifshitz 
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argument, applicable to the low energy DoS for a low concentration of scalar impurity, leads 
instead to 39 a behaviour 

PLif.(^)~e- C2 ^, (114) 

where C2 is a constant. Our choice for this fit is motivated by a recent numerical work |184j 
concerning the area A of the outer boundary of planar random loops. In this work, the author 
suggests that the limit distribution of A is the Airy distribution implying, for small A values, 
a behaviour of the type exp(— const. /^4 2 )/^4 2 . Remarking that A = J2 n even + Yin odd^™' 
it is natural to expect for the distribution of the random variable S, eq. (jlllj) . a behaviour at 
small S that is roughly given by exp(— const. /S 2 ). Thus, we deduce the form pat(E) for the 
low-energy fit of (p(E)). 

• Transition between a. = 0.5 and a — ► 0. 

Finally, when a grows from to 0.5, the oscillations in the spectrum must disappear at some 
critical value a c (see figure ITT]) . Numerical simulations, specific heat considerations and, also, 
diagrammatic expansions [77] give a c ~ 0.3. 




012345678 012345678 
EI<B> EI<B> 

Figure 11: Average density of states of the Hamiltonian U00\) for different values of the flux per 
tube 4> = oi(f)Q. The average magnetic field reads (B) = 2ir[ia. We recall that the corresponding 
Landau spectrum is given by E n = (B) (n + 1), for n 0. From ref. 



39 The Lifshitz argument applies to the DoS for the random Hamiltonian H = — A + X^i u{r — ri) in dimension d, 
where u(r) is a sharply peaked scalar potential. For a low density of impurities fi, the DoS behaves as Plh.(E) ~ 
exp(— a fi/E dl12 ) at low energy. 
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7 Conclusion 



A non experienced reader may have the feeling that the two topics which have been covered in this 
review, namely one-dimensional disordered systems and quantum graphs, are essentially disjoint. 
In fact there are many interesting links between these two topics, both at a methodological and 
a conceptual level. The use of metric graphs for modelling quantum phenomena observed in 
disordered metals goes back to the pioneering work of Shapiro |191j and Chalker & Coddington 
|56j . These systems have been used to study quantum localization and more recently as a model 
system for spectral statistics (for a recent review, see ref. [202 J. Another set of similar questions 
is provided by the study of scattering properties of chaotic graphs |143| and disordered systems 
|1()3| (see also section H~2"l and ref. |174j for recent developments). We have also seen that apart 
from its interest to study the spectral properties of metric graphs, the spectral determinant 
also allows studying several properties of networks of quasi-one-dimensional weakly disordered 
wires [T771 IT7SI 13] . 

As a conclusion we would like to mention several open problems : 

• A study of spectral statistics in the case of graphs with a random Schrodinger operator [there 
is still a factorised structure but the matrix R is now given by eq. (|66j)] . 

• A probabilistic understanding of the star graphs using the tools of excursion theory developed 
by Barlow et al. 23 in the context of the Brownian spider. A first step would be to recover 
those probabilistic results (the joint law of the occupation time inside the branches) by a spectral 
approach. It is however not excluded that the probabilistic approach could provide a key to 
a deeper understanding of those quantum systems. In the context of classical systems such 
probabilistic approaches have been very useful, e.g. recent studies of the Stochastic Loewner 
Equation have made enormous progress in understanding the statistical physics of a class of 
two-dimensional systems. This calls for new probabilistic techniques for quantum systems as 
well. 

• Several functionals of the Brownian motion (|13I16() appear when studying the important 
question of dephasing due to electron-electron interaction in networks of quasi-one-dimensional 
weakly disordered wires. The fact that such simple functionals appear relies on the translation 
invariance of the two particular problems studied (an infinite wire (2] and an isolated ring |212j ). 
For a network with arbitrary topology, the relevant functional of the Brownian bridge x(t) 
is given by f£ dr W(x(t), x(t — r)) where W(x,x') = ^[Pd(x,x) + Pd(x',x')] — Pd(x,x') with 
—APd(x,x') = 5{x — x'). It now involves a nonlocal functional in time which is difficult to 
handle. Progress in this direction would allow clarifying the interplay between the electron- 
electron interaction and the geometrical effect and help in analyzing recent experimental results 

inn uni 

• The question of extreme value spectral statistics was addressed in the framework of random 
matrix theory |214| and these studies have found several applications in the context of out-of- 
equilibrium statistical physics (see the review |215j ) . However this question was first addressed in 
the context of one-dimensional disordered systems |117j . The study of supersymmetric random 
Hamiltonian, for which the bottom of the spectrum plays a special role, has emphasized the 
interest in extreme value spectral statistics |2()4| (in particular it indicates level correlations). It 
would be interesting to extend such studies to other models and find some physical situations 
where these results would be applicable. 

• A beautiful heuristic argument was provided by Lifshitz to explain the spectral singularity 
for Hamiltonians with a weak concentration of localized scalar impurities (see for example the 
book |157| 1. Other mechanisms should be invoked to explain the nature of the quantum states 
responsible for the low energy power law behaviour in the case of a uniform strong magnetic 
field with 5-impurities jlUUI HUf] . For the model of randomly distributed magnetic fluxes, the 
numerical simulations for a = 1/2 have suggested the new type of singular behaviour (|113j) . It 
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would be interesting to provide heuristic arguments to understand more deeply the origin of this 
singularity. 
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